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Abstract. Since the beginning of the past century, the nature of HD 108 has been a subject of intense debate. One after another, 
astronomers explored its variability and attributed it either to binarity, or to changes in the stellar wind of a single star In this 
article, we analyse a 30 years long campaign of spectroscopic observations of this star with special emphasis on the last 15 
years during which photographic plates have been replaced by CCD detectors. Our investigation of the radial velocities of 
HD 108 yields no significant short- or long-term period and does not confirm the published periodicities either Though the 
radial velocity of HD 108 appears clearly variable, the variations cannot be explained by the orbital motion in a spectroscopic 
binary. However, our data reveal spectacular changes in the H I Balmer lines and some He I profiles over the years. These lines 
continuously evolved from P Cygni profiles to 'pure' absorption lines. A similar behaviour was already observed in the past, 
suggesting that these changes are recurrent. HD 108 seems to share several characteristics of Oe stars and we discuss different 
hypotheses for the origin of the observed long-term variations. As we are now in a transition period, a continuous monitoring 
of HD 108 should be considered for the next few years. 
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1. Introduction 

HD 108 has stimulated a lot of controversy over the past 75 
years. One of the first detailed studies of its spectrum was 
presented by Andrillat et al. (1973), who also reviewed the line 
profile changes (from P Cygni to pure emission/absorption 
lines). Some years later, using the radial velocities of the He II 
AA4200, 4542 absorption lines, Hutchings (1975) derived an 
orbital period of 4.612 days, but this result was not confirmed 
by Vreux & Conti ( |1979[ ). Vreux & Conti ( |1979| ) suggested 
that the radial velocity variations were rather due to wind 
variability and though they pointed out the possibility of a 1.02 
day period, they concluded that the variability of HD 108 was 
most probab ly not periodic. On the basis of Hutchings's work, 
Bekenstein ( 1976 ) suggested that HD 108 had undergone a 
type II supernova without breakup of the binary. Combining 
radial velocity data from earlier st udies with their own mea- 
surements, Aslanov & Barannikov (1989 1 argued in favor of a 
5.7937 days period, again attributed to binarity. But the debate 
about the re al nat ure of HD 108 was not settled. More recently. 
Underbill (1994) found no periodic variation of the radial 
velocities, thus rejecting the binary hypothesis. She suggested 
instead that HD 108 was a single star surrounded by a disc and 



displaying jets that were alm ost per pendicular to the line of 
sight. According to Underbill (1994), this particular geometry 
accounts for the peculiar line profile morphology and could 
explain the spectral variations seen in HD 108. Using V-band 
polarimetric observations. Fox & Hines ( 1998| ) did not find 
any evidence of a binary nature of HD 108, but their results 
did not allow them to choose between large-scale jets or 
small-scale blobs as the structures triggering the variability. 



On the contrary, Barannikov (1999) combining his own recent 



radial velocity measurements with all older ones, claimed the 
existence of a longer period of 1627.6 days and concluded 
again that the star is indeed a binary. 

Other properties of HD 108 are also intriguing : its UV spec- 
trum suggests a low mass-loss rate (Hutchings & van Heteren 



1981), incompatible with the one derived from its optical and 



IR (Ferrari-Toniolo et al. 1981) spectrum. Moreover, some 
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authors marked HD 108 as a runaway star, though they did not 
agree on the exact value of its peculiar velocity. 

In this paper, we will analyse CCD spectra collected over 
fifteen years. The observations are presented in Section 2, and 
we investigate the main spectral characteristics of HD 108 in 
Section 3. In Section 4, we will discuss the short- and long- 
term variations of the radial velocities and search for a periodic 
behaviour of HD 108. In the next sections, we analyse the vari- 
ability of the equivalent widths and the line profiles. Finally, we 
examine the photometric data and the proper motion of HD 108 
before we conclude in Section 8. 
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2. Observations and Data Reduction 



Table 1. Summary of our CCD spectroscopic observations. 
indicates the total number of spectra obtained during a run. 



Date 


Telesc. 


Wav. range 


N 


S/N 


Disp. 












(Amm^^) 


July 1986 


1.93m 


3875-4320 


2 


275 


33 


July 1987 


1.93m 


3940-4400 


20 


280 


33 


Aug. 1987 


1.93m 


4415-4915 


I 


200 


33 


Aug. 1989 


1.93m 


3875-4915 


I 


300 


130 


Sep. 1989 


1.93m 


8430-10400 


I 


25 


260 




1.93m 


8440-11130 


I 


290 


260 


Aug. 1990 


1.93m 


8450-8765 


1 


240 


33 




1.93m 


6830-7130 




150 


33 




1.93m 


6840-9945 


2 


600 


260 


Aug. 1991 


1.93m 


4250-4660 


6 


360 


33 


Oct. 1992 


1.93m 


8450-8765 


3 


200 


33 


Oct. 1993 


1.93m 


3940-4360 


2 


220 


33 


Aug. 1994 


1.93m 


3930-4360 


7 


530 


33 


Aug. 1996 


1.52m 


4065-4925 


3 


230 


33 


Feb. 1997 


1.52m 


6500-6710 


4 


150 


8 


July 1997 


1.52m 


4100-4960 


3 


280 


33 


Sep. 1998 


1.52m 


4455-4885 


3 


400 


16 


Nov. 1998 


1.52m 


4435-4560 


15 


100 


5 


July 1999 


1.52m 


4060-4930 


3 


300 


33 


Aug. 1999 


1.52m 


4060-4930 


4 


300 


33 


Sep. 2000 


1.52m 


4455-4905 


11 


300 


16 



An extensive spectroscopic survey of HD 108 was con- 
ducted from 1986 to 2000 with the Carelec and the Aurelie 
spectrographs, fed respectively by the 1.93 m and 1.52 m tele- 
scopes of the Observatoire de Haute-Provence. For Carelec, 
the detector used was a thin back illuminated RCA CCD 
with 323x512 pixels of 30/im squared. For Aurelie, it was a 
TH7832 linear array with a pixel size of 13/im until 1999. In 
2000, this detector was replaced by a 2048 x 1024 CCD EEV 
42-20#3, whose pixel size is 13.5/im squared. The exact wave- 
length ranges and mean S/N are given in Table |l], together with 
the dispersion. All the data were reduced in the standard way 
using the IHAP and MIDAS softwares developed at ESO. The 
spectra were normalized by fitting splines through carefully 
chosen continuum windows. Some older photographic plates 
taken in 1971-72, 1974 and 1976 at OHP were also used to 
study the line profile variability. 

3. Tlie spectrum of IHD 108 

The spectrum of HD 108 is quite peculiar, showing a great 
number of emission lines. The blue part of the spectrum is 
dominated by the presence of strong N III A A 4634-41 and 
He II A 4686 emission lines (see Fig. |l]). As for many Of stars 
(Underhill et al. 1989), these lines stand on top of a broad 
emission bump between 4600 and 4700A. In addition. Si III 
AA4552, 4568, 4574 and various O II emission lines, as well 
as the unidentified Of emissions situated near AA4485, 4504 
are also present. All Balmer lines except Ha appear as P Cygni 
profiles during most of our observing runs, but their morphol- 



ogy changes with time (see Section 6). On the contrary, some 
lines always appear in absorption, for example Si IV AA4088, 
4116 and He II AA4200, 4542. Some interstellar features can 
also be seen in Fig |l]: besides the Ca II H and K lines, there are 
diffuse interstellar bands (DIB's) around 4428, 4726 and 4763 
A, as well as a CH Une at 4300 A. 

The Ha and He I AA6678, 7065 emissions dominate the 
red part of the spectrum. In our 1989 near-IR spectrum, the 
Paschen lines of Hydrogen appear in emission, as well as He I 
A 10830 which is very intense. He II A 10124 is also present, 
but in absorption. 



3.1. Spectral type 

HD 108 has been classified successively as 06fqkQ (Plaskett 
& Pearce |93T|), 08fp (Morgan et al. [19551 ), 07If (Conti & 
Alschuler 1971 ), 06:f?pe (Walborn |1972[) and finally 07 Illfpe 



(Underhill 1994) 



Applying Conti & Alschuler's (1971 ) criterion to our data, the 



spectral type varies between 04 (Aug. 1987) and 07.5 (Sept. 
2000). A similar type of variations has already been reported 



by Reals ( |1950[ ), who found for HD 108 an 06 type in 1938 
and 07 in 1945. A closer inspection of this situation indicates 
that the He II A 4542 line profile is quite constant while the 
spectral type variations are due to the varying shape of the He I 
A 4471 line which sometimes appears with a P Cygni profile 
(see Sections 5 & 6), revealing that this line is not completely 
of photospheric origin. This casts doubt on the possibility to 
apply the Conti & Alschuler criterion to this star. 
For the luminosity class, we use the criterion based on the value 
of \og(W") = log wfHcTA^4i°43)^ ' The equivalent widths 
(EWs) of these two absorption lines show no strong variations, 
and result in a supergiant classification, regardless of the year 
of observation. Even though it is possible that He I A 4143 
could also be contaminated by a variable emission component, 
this effect should be small with regard to the stabili ty of our 
result. Using polarimetric data. Fox & Hines (1998) also fa- 
vored a supergiant classification. The presence of a strong He I 
A 10830 emission is also consistent with a supergiant classifica- 
tion, although this emission is also observed in Oe-type objects 
(Vreux & AndrilIat[T979|). 



The presence of N III A A 4634-41 and He II A 4686 in emission 
justifies an 'f tag. Moreover, a 'p' tag can be added because 
of the presence of many emission lines, more numerous than 
those commonly found in Of spectra. Finally, the emission in 
H I Balmer lines results in the addition of an 'e' tag. Choosing 
the latest spectra, i.e. the least affected by possible emission 
in He I, we can derive a probable 07.5 Ifpe type for HD 108. 
If the He I A 4471 absorption continues to strengthen over the 
next years (see Sections 5 & 6), it is possible that we will finally 



recover an 08 type as found by Morgan et al. (1955) 



' The T tag stands for N III AA4634, 4641 and He II A 4686 
present in emission, the 'q' is for the presence of P Cygni profiles 
and 'k' for the presence of Ca lines which are of interstellar nature. 
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Fig. 1. Composite spectrum of HD 108 from blue to IR : the AA 3875-3942 domain comes from the 1986 data, 3942-4300 from 
the 1994 data, 4300-4930 from the 1999 data, 6500-6710 from the 1997 data, 7000-7100 from the 1990 data and 8430-11000 
from the 1989 data. 



4. Is HD 108 a spectroscopic binary? 

Using the rest wavelengths listed in Table ||, we determined 
the radial velocities by fitting gaussians to the bottom (top) end 
of the absorption (emission) lines. The radial velocities (RVs) 
of the He II A A 4200, 4542 lines are summarized in Table ^ 
together with the RVs of the N III A A 4634,4641 and He II 
A 4686 emission lines. We hav e used the generalized spectro- 
gram technique of Heck et al. (1985, hereafter HMM) and the 
trial method of Lafler & Kinman ( 1965 ) to search for periodic- 



ities in the He II AA4200, 4542 RV time series. For the HMM 
method, we evaluate a rough estimate of the significance level 
(SL) by SL = e ^ where P is the periodogram amplitude and 
cr the rms dispersion of the data. We will consider that a period 
is significant if its SL is less than 0.01. 

4. 1. Short-term variability 

Over all those years of intense study of HD 108, only a few au- 
thors really addressed the issue of short-term variability of this 
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Table 2. Mean heliocentric radial velocities of He II AA4200, 4542. We have also included the RVs of N III AA4634, 4641 and 
He II A 4686 for completeness. The date of the observation is given in the format HJD-2440000. 



Run 


Date 


RV4200/4542 


RV4634/4641 


RV4686 


Run 


Date 


RV4200/4542 


RV4634/4641 


RV4686 






km s^^ 


km 


km 






kms~^ 


km 


km 


Jul. 1986 


6627.582 


-66.4 






Aug. 1996 


10318.644 


-56.8 


-66.8 


-53.9 




6628.548 


-47.0 






Jul. 1997 


10638.574 


-60.8 


-63.0 


-41.4 


Jul. 1987 


7007.536 


-68.5 








10639.557 


-61.7 


-61.5 


-54.8 




7007.544 


-79.2 








10642.573 


-59.1 


-60.3 


-53.4 




7007.549 


-72.9 






Sep. 1998 


11066.635 


-64.1 


-77.0 


-65.6 




7007.554 


-74.2 








11070.629 


-65.0 


-76.8 


-64.9 




7008.502 


-84.4 








11072.638 


-62.8 


-76.6 


-61.7 




7008.517 


-80.4 






Nov. 1998 


11132.377 


-73.1 








7008.524 


-84.5 








11132.421 


-66.5 








7009.526 


-82.9 








11133.330 


-69.4 








7009.538 


-78.4 








11133.380 


-68.2 








7010.549 


-74.8 








11133.432 


-69.7 








7010.558 


-73.4 








11133.589 


-71.6 








7011.554 


-72.8 








11134.378 


-67.9 








7011.561 


-75.8 








11134.448 


-66.9 








7011.568 


-77.4 








11135.339 


-68.4 








7011.572 


-74.9 








11135.394 


-69.8 








7012.566 


-77.2 








11135.458 


-66.4 








7012.573 


-71.9 








11136.373 


-67.0 








7013.540 


-70.7 








11136.439 


-69.5 








7013.545 


-64.7 








11137.401 


-67.6 








7013.549 


-70.2 








11137.467 


-67.8 






Aug. 1987 


7016.601 


-69.2 


-75.4 


-72.2 


Jul. 1999 


11374.597 


-70.1 


-75.7 


-68.5 


Aug. 1991 


8490.478 


-93.5 


-93.3 






11376.532 


-75.7 


-78.2 


-69.6 




8490.489 


-86.4 


-85.1 






11378.530 


-71.2 


-75.8 


-65.1 




8490.505 


-85.2 


-84.0 




Aug. 1999 


11396.600 


-68.8 


-74.4 


-66.4 




8492.739 


-94.2 


-89.9 






11403.635 


-73.4 


-79.2 


-68.9 




8496.590 


-84.6 


-80.4 






11406.634 


-69.7 


-76.0 


-67.1 




8497.581 


-77.0 


-72.6 






11407.639 


-71.6 


-77.3 


-69.4 


Oct. 1993 


9264.447 


-61.4 






Sep. 2000 


11810.586 


-69.4 


-79.6 


-66.9 




9264.452 


-67.3 








11810.595 


-70.6 


-79.6 


-67.5 


Aug. 1994 


9576.553 


-44.2 








11811.567 


-68.6 


-78.5 


-67.4 




9577.542 


-64.0 








11811.575 


-68.6 


-78.9 


-67.2 




9579.522 


-64.8 








11812.594 


-71.6 


-80.0 


-70.4 




9580.481 


-58.0 








11813.603 


-67.8 


-76.5 


-67.8 




9581.535 


-58.8 








11814.595 


-69.7 


-79.4 


-66.2 




9582.511 


-63.8 








11815.596 


-69.8 


-79.8 


-69.4 




9583.533 


-59.3 








11818.531 


-66.3 


-77.9 


-64.4 


Aug. 1996 


10316.650 


-58.9 


-63.4 


-56.9 




11819.560 


-68.3 


-79.2 


-66.7 




10316.658 


-58.3 


-63.0 


-57.0 




11821.561 


-72.0 


-80.7 


-69.3 



Table 3. Adopted rest wavelengths of the lines used to compute 
the RVs displayed in Figs.|^, || and ^ and listed in Table I 



Ion 


Effective 


Nature 




wavelength 




He II 


4199.830 


Abs. 


Hell 


4541.590 


Abs. 


Si III 


4552.654 


Em. 


Si III 


4567.872 


Em. 


NlII 


4634.250 


Em. 


NlII 


4641.020 


Em. 


Hell 


4685.682 


Em. 



star, even if most of the periods found were rather short. For 



instance, Vreux & Conti ( 1979 ) obtained two spectra per night 
during two consecutive days of their observing run. However, 
to analyse completely the short-term variability, more spectra 
taken on numerous consecutive days are needed. To investigate 
this problem, we collected at least two spectra per night during 
the 7 days of the August 1987 run and during the 6 days of the 
November 1998 run. The 1987 data cover the He II A 4200 line, 
while the 1998 data cover the He II A 4542 line. These data en- 
able us to search for short-term variations of the spectra, and 
also to check the short periods (a few days) given in the liter- 
ature. Because of the poorer spectral resolution (1.7 A) of the 
1987 data, the RVs show a large scatter in 1987 compared to 
the 1998 measurements that were obtained with a resolution of 
O.2A (see Figs. I and |). 
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IIuLchmgs's period : 4.6117 d 

Fig. 2. RVs of He II A 4200 as measured on our 1987 spectra 
(open squares) and RVs of He II A 4542 from our 1998 data 
(filled triangles) folded with Hutchings' period of 4.61 17 d and 
superimposed on his orbital solution. 



We tried to search for a periodicity in the He II RVs mea- 
sured on the two sets of data (20 spectra in 7 days +15 spec- 
tra in 6 days), but we did not find any significant short pe- 
riod. Moreover, those data do not confirm the periods given 
in the literature. In Figures ^ and|^, our data are folded respec- 
tively with Hutchings'period (4.6117 d, Hutchings 1975) and 
the 5.7937 days period of Aslanov & Barannikov (1989). Both 
periods can clearly be ruled out. Due to the severe aliasing of 
our time series, the 1 .02 day period taken from Vreux & Conti 
(1979) cannot be sampled over an entire cycle with our data. 
Though the scatter of the folded data (20 km s^^ for A(p ^0. 1 
in 1987 and 10 km s^^ for Acf) ~0.3 in 1998) argues against 
this last period, it cannot be totally excluded. Meanwhile, we 
emphasize that no short period, typically from 1 to 5 days, does 
appear in our data. 

We have also applied the period research techniques to the He II 
EWs determined in those two runs. Again, no significant period 
appears. With the 1998 high-resolution data, we also searched 
for short-term variations in the line profiles: neither He I A 447 1 
nor He II A 4542 seem to change during the whole duration of 
our run. Short term variations of typically 1 to 5 days - if they 
exist - are thus only of very small amplitude. 



-100 



4.2. Long-term behaviour 

Since we do not detect any short-term variability in the He II 
RVs, we can average the He II absorption RVs for each observ- 
ing run: we have plotted those means in Fig. and we have 
combined these values with all published RVs in Fig. ^. We 
caution that the oldest data are rather uncertain, and that they 
should only be used for a qualitative analysis. 



Aslanov & Baranriikov's period : 5.7937 d 

Fig. 3. Same as Fig. |[ but for the period of 5.7937 d proposed 
by Aslanov & Barannikov (1989). 



The He II AA 4200, 4542 lines are most probably formed in 
the deeper (photospheric) layers of the stellar atmosphere and 
can thus give rather unpolluted information about the star's ac- 
tual motion. Other absorption lines (Mg II A 4481, Si IV A 4631 
and all N III lines between 4510 and 4534 A) and some emis- 
sion lines (the unidentified Of emissions AA 4486, 4504 ; Si III 
A A 4552, 4568 ; N iii A A 4634, 4341 ; He ii A 4686 and O ii 
A 4705) display the same behaviour as the He II absorption 
lines. The RVs of some of these lines are shown in Fig. ^ On 
the other hand, the RVs of the main absorption components of 
H/3, H7 and He I AA4388, 4471, 4713 become less and less 
negative over all the years of our observing campaign. This 
is linked to the long-term line profile variations discussed in 
Section 6. 

The continuous trend of the absorption lines towards mo re neg - 
ative RVs since 1922-1923, as reported by Underhil^ ( |l994| ). 
does not appear in our data. On the contrary, in 1994 and in 
1996, the RV of He II A A 4200, 4542 is about the same as in 
1922. The absence of a clear very long term trend is best seen 
in Fig. ^ 

Underbill ( |1994 ) also mentionned that the RVs of the absorp- 
tion and emission lines show different behaviours. But in our 
data. He II absorptions and He II, O II, Si IV and N III emis- 
sions display quite similar behaviours (as can be seen on Fig. 
^ and on Fig. 1 of Barannikov 1999[ ). However, the 'abrupt' 
shift of the RVs towards more negative values in 1991, seen by 
Underbill (1994) does appear in our data as well as in those of 
Barannikov (1999). If intrinsic to a single star, this behaviour 
could indicate an episodic change of the velocity structure in 



" The errors quoted in Table 3 of Underbill ( 1994 ) are not tbe stan- 
dard deviations of tbe data, contrary to what we used in Fig § or ^ 
but tbe standard deviations of tbe means: those two values differ by a 
factor \//V. This explains the rather low value of the a in ber Table. 
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Fig. 4. Mean RVs measured for each observing run: there are 
He II AA 4200, 4542 absorptions (represented by filled trian- 
gles) and some emission lines like Si III AA 4552, 4568 (stars), 
N III AA4634, 4641 (hexagons) and He II A 4686 (asterisks). 
The bottom axis yields the date in heliocentric Julian days 
(HJD), while the x-axis on top of the panel indicates the year of 
observation: tickmarks represent the middle of the year quoted 
above. 1-cr error bars are shown for He II AA4200, 4542 RVs 
(when no error estimates are shown, it means that only one 
spectrum was available). 



the inner regions of the expanding atmosphere. 
Using our data set only, a period of about 4600 days appears in 
the Fourier periodogram. However, with 5i — 0.89, this detec- 
tion is not significant. Moreover, this period completely disap- 
pears when we include the RVs determined by Plaskett (1924 ), 



Hutchings ( 1 975 ), Vreux & Conti ( 1 979 1, Underbill ( 1 994 ) and 



Barannikov (1999). Considering the whole data set, no signifi- 



cant period becomes visible, the minimum SL reaching 0.93. 



According to Barannikov (1999), the He II RVs should follow 



a period of 1627.6 d. But when we fold our data according to 
this period, it clearly appears that it is ruled out by our data (see 
Fig. Barannikov used low-resolution^ photographic spec- 
tra, which are prone to errors (there are differences up to 70 
kms^^ between two consecutive nights for a 1627.6 days pe- 
riod with an amplitude K of only 10.5 kms^^!): this might 
explain the important scatter seen in his orbital solution. We 
note that the minimum in Barannikov's radial velocity curve 
occurred in 1991, i.e. at the same epoch as we observe the RV 
shift discussed hereabove. However, our data do not suggest 
that this phenomenon is periodic. 

5. Equivalent Widths 



Aslanov & Barannikov (1989) claim that the equivalent width 
variability of emission lines is of stochastic nature, while the 

^ The dispersion was 44 A/mm^'^. 
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Fig. 5. All published radial velocities of absorption lines 
(mainly He II A 4200, 4542): the open squares represent 



the mean data from Plaskett (1924), the filled squares from 



Hutchings (1975), the open triangles from Vreux & Conti 



(1979), the stai-s from Underbill (1994) and the filled triangles 



from this work. 



one associated with absorption lines exhibits a period twice 
the one derived from their radial velocities. More recently. 
Underbill (1994) found no periodicity in the variability of the 
EWs. 

In our data, two different categories of lines emerge: most of 
the lines appear rather constant but several others have clearly 
variable EWs. In the first category. Si IV A 4088, He II A 4200, 
4542 and He I A 4143 exhibit maximum EW variations below 
^ 20% and the scatter is quite low (see Table Other lines fall 
also into this category, e.g. the O II and C II emission lines, the 
unidentified Of emissions and the N III absorptions. In agree- 
ment with Underbill ( 1994 ), we do not find any significant peri- 
odicity when considering the variability of the EWs of the He II 
absorption lines. 

Turning now to the lines showing variable EWs, we will fo- 
cus here on the behaviour of those lines for which we have the 
best coverage with our data sets. The Si III AA4552, 4568 and 
C III AA 4647-50 emissions seem to decrease gradually during 
the whole duration of our observing campaign. The EWs are 
reduced by a factor 2 between 1 987 and 2000. The He I A 47 1 3 
line also presents EW variations, up to 50% from the mean. 
But the most impressive changes occur in H7, He I A 4471 and 
H/3. Their total equivalent widths (absorption H- emission) in- 
crease gradually (see Table In Fig. ^ one can clearly see 
this continuous trend to higher EWs, as the emission compo- 
nent weakens and finally disappears (see Section 6). 
This separation of the lines in 2 categories, the rather constant 
ones (e.g. He II emission and absorptions) and the clearly vari- 
able ones (e.g. Si III and the most variable He I and H I lines). 
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1988 1990 1992 



99/1 1996 1998 2000 



Fig. 6. Our mean He II data folded with Barannikov's (1999) 
period of 1627.6 d, superimposed on his orbital solution. 1-a 
error bars are shown: when no error estimates are shown, it 
means that only one spectrum was available. 



points most probably towards different formation regions for 
each category. 

Table 4. Mean EW for different lines, averaged over the entire 
observing campaign. The quoted EWs for H/3 and H7 are the 
total ones, including the He II Pickering absorption (and the 
DIB in the case of H/3). The last two columns summarize the 
results of a TVS analysis (see Section 6). The VTVS is given 
in units of the continuum flux, while the last column provides 
the wavelength interval over which the detected variations are 
significant at the 99% level. 



Line 


EW(A) 


VTVS 


width(TVS) 




(A) 




(A) 


Si IV A4088 


0.196±0.021 






He I A4143 


0.040±0.006 






He II A4200 


0.538±0.033 






H7 


0.506±0.385 






He I A4471 


0.643±0.215 


0.065 


4462.3-4473.3 


He II A4542 


0.672±0.057 


0.017 


4535.5-4542.6 


Si III A4552 


-0.141± 0.045 


0.016 


4549.9-4555.3 


Si III A4568 


-O.llOi 0.030 


0.010 


4565.4-4569.4 


N III A4634 


-1.040±0.059 


0.014 


4630.7-4637.0 


N III A4641 




0.018 


4637.0-4644.4 


C III AA4647-50 


-0.306±0.077 


0.013 


4644.4-4652.7 


He II A4686 


-0.524±0.045 


0.013 


4709.2-4714.2 


Hel A4713 


0.134±0.040 


0.092 


4849.0-4866.5 


H/3 


1.457±0.788 








7000 8000 9000 lo* 

HJD - 2 440 000 



1.1x10' 



i.axio" 



Fig. 7. Average EW for each observing run: He II A 4200 is 
represented by filled triangles, H7 by open squares. He I A 447 1 
by crosses and H/3 by open triangles. The errors are very small, 
roughly the size of the symbols used in the figure, and they are 
listed in Table || 



6. Line Profile Variabiiity 



As summarized for example in Andrillat et al. (1973), many 
line profiles of the spectrum of HD 108 vary over the years, 
with some of them changing spectacularly from P Cygni pro- 
files to pure absorption or emission lines. This is especially 
the case for the Balmer Hydrogen lines and some He I lines 
(AA4388, 4471, 4713). In our data. He l AA4388, 4713 ap- 
peared only once, in 1987, with a P Cygni profile while these 
lines display pure absorptions since then. 
To better quantify the line variability, we used the temporal 
variance spectrum (TVS), as defined by Fullerton et al. ( 1996 ). 
This TVS analysis, when applied on the common spectral range 
of the data collected in 1987 and from 1996 to 2000, yields in- 
teresting results (see Table ^: in addition to the He I and H I 
lines, the absorption fines of N III (from 4510A to 4534A) and 
of He II A 4542, as well as the emission fines of N III AA 4634- 
41, He II A4686, Si III AA4552, 4568 and C III AA 4647-50, 
do all vary, though with a smaller amplitude. On the contrary, 
other emission lines (e.g. O II and the unidentified AA4486, 
4504) were stable during the whole observing campaign. 
To compare the variations of the different lines, we have 
used the local pattern cross-correlation technique discussed by 
Vreux et al. ( |1992| ) and by Gosset et al. ( |1994j ). But before that, 
we have convolved the 1998 and 2000 data until they reach 
the same resolution as the 1996, 1997 and 1999 spectra. Only 
then, we have cross-correlated the variation pattern of some un- 
blended lines (e.g. He I A 4471, He II A 4686) with the whole 
spectra. Using the He II A 4686 emission line as a reference, 
a strong correlation appears between the deformation pattern 
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Table 5. Average EW for each observing run for He II A 4200, H7, He I A 4471 and H/3. Standard deviations are also quoted : 
when no error estimates are given, it means that only one spectrum was available. 



Run 




EW He II 




EWH7 


EWHel 




EWH/3 


Jul. 1986 





559 ± 0.003 












Jul. 1987 





524 ±0.019 





153 ± 0.038 








Aug. 1987 










0.120 




0.566 


Aug. 1991 









367 ± 0.034 


0.235 ±0.019 






Oct. 1993 





548 ± 0.051 





577 ± 0.075 








Aug. 1994 





518 ±0.012 





635 ±0.017 








Aug. 1996 





557 ± 0.049 





781 ±0.017 


0.417 ±0.017 





610 ±0.066 


Jul. 1997 





562 ± 0.033 





830 ± 0.045 


0.542 ± 0.006 





745 ± 0.018 


Sep. 1998 










0.622 ± 0.029 






Nov. 1998 










0.696 ± 0.039 






Jul. 1999 





539 ±0.018 


1 


221 ± 0.055 


0.694 ±0.021 


1 


267 ± 0.070 


Aug. 1999 





587 ± 0.057 


1 


232 ± 0.043 


0.743 ±0.017 


1 


155 ± 0.045 


Sep. 2000 










0.881 ± 0.036 


2 


229 ± 0.090 


Mean 


0.538 ±0.033 


0.506 ± 0.385 


0.643 ±0.215 


1.457 ± 0.788 



of He II A 4686 and the ones of the N III AA4634, 4641 Hnes. 
On the contrary, the He I A 4471 pattern seems rather related to 
the variation of the N III absorptions between 45 10 and 4534A, 
He II A 4542, C III A 4647-50, He l A 47 13 and H/5. The be- 
haviour of the N III A A 4634-41 lines suggests therefore that 
these lines form in the same physical region as the He II A 4686 
emission, while the clearly distinct variations of the He I and 
Balmer emissions point towards a different origin and/or emis- 
sion mechanism for the latter lines. The variability pattern of 
He II A 4686 and N III A A 4634-41 between 1996 and 2000 re- 
veals an enhancement of the core of the line while the red wing 
is progressively depleted. 

These N III profile variations probably point towards a wind 
contribution to the emission. In this context, it is worth recall- 
ing that once that a velocity gradient exists, the fluorescence 
mechanism, suggested by Swings ( |1948t ) to account for the oc- 
currence of the N III A A 4634-41 emission lines in the spectra 
of Of stars, can become extremely effective (Mihalas 1973). 
The most impressive variability is unquestionably displayed by 
Ha, H/3, H7 and He I A 4471, as already quoted in Section 5. 
The long-term line profile variability of these lines is presented 
in Tableland Fig.^ the changes are clearly outstanding. We 
do not have a lot of information about Ha, but the line is clearly 
variable. Our 1997 Ha data yield a normalized peak intensity 



2 - 



of 1.8, whereas in 1990 (Underhill 1994), it reached a value 
of 2.5. A similar behaviour was observed by Beals (1950): in 
1938, the normalized intensity of Ha was over 2, but it was 
only 1 .3 in 1945. The Ha line thus roughly behaves in the same 
way as H/3, H7 and He I A 4471. Other Balmer lines (e.g. HS, 
He) and He I lines (e.g. 4388, 4713) also display a similar be- 
haviour, but we will now focus on the H7, He I A 4471 and H/3 
lines, for which we have the most extensive data set. 
These three lines displayed a P Cygni profile at the beginning 
of our observations in 1987. Then the emission component pro- 
gressively weakened until the lines appeared completely in ab- 
sorption. The first line that went into absorption is He I A 4471, 
followed by H7 and finally H/3. After this change, the depth 
of the He I absorption lines continued to increase gradually, as 



1.6 



1 1 1 1 1 1 1 1 


\1987 




V 1995 




V^J998 




^1999 




I 2000 





4620 



4640 

Nlll. CI 



&■ lie 



4660 

I omission 



Fig. 8. Aspect of N III AA 4634-41, C III AA 4647-50 and He II 
A 4686 for different years: the dispersion used in 1987, 1996 
and 1999 is 33 A mm-\ while it is 16 A mm"^ in 1998 and 
2000 (see Table |l]). The C III lines show a very significant de- 
crease in intensity between 1987 and the recent years, while all 
lines seem sharper in 2000. The heliocentric correction has not 
been applied to the spectra displayed on this figure. 

can be seen in Figs. ^ and ^. Since the change into absorption 
occurred only recently, we have no information yet about the 
behaviour of the H7 and H/3 lines after the emission has disap- 
peared. 

Interestingly, a similar transition was observed some 56 years 
ago (see Table ^: all three lines were seen as pure absorptions 
in 1944 and also in 2000. In 1934, He I A 4471 was already 
described as an absorption, while H7 and H/3 were still dis- 
playing P Cygni profiles. Then, H7, followed very soon by H/3, 
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went into absorption. Finally, in 1954, a weak emission reap- 
peared in H/3. Therefore, the long term line profile variations 
of HD 108 seem to be recurrent with a timescale of about 56 
years. 

However, the interpretation of the Balmer line profile variations 
is complicated by the fact that these lines are blended with He II 
Pickering lines. Following the same procedure as in Lamers & 
Leitherer (1993), we attempted to restore the H I line profiles 
by subtracting a fake He II line obtained by interpolation (H7) 
and extrapolation (H/3) from the unblended He II A 4200 and 
A 4542 lines. The reconstructed profiles are striking: the hydro- 
gen lines no longer appear as P Cygni profiles, but as rather 
sharp and symmetrical emissions superimposed on broader ab- 
sorptions (see the lower panels of Fig. This emission fades 
with time, and finally nearly disappears in H7 in 1999. The 
RVs of the 'restored' H/3 and H7 emission peaks become more 
negative during the last two years. 

We have also tried to remove the photospheric He l absorp- 
tion to recover the stellar wind line profile using two slightly 
different approaches. First, we assume that the emission has 
completely disappeared in the spectra observed in 2000, and 
use this profile for the He I photospheric contribution. Second, 
we evaluate the photospheric absorption by fitting a gaussian 
profile to the blue wing of the A 4471 line, assuming it is not 
contaminated by emission. In this latter case, we set the RV of 
the photospheric He I component to be equal to the observed 
RV of the He II AA 4200, 4542 lines. The resulting pure emis- 
sion profiles are displayed in Fig. ^. We caution however that 
the reconstruction of the He I A 447 1 wind profile is far more 
uncertain than the restoration of Balmer lines. 
The rather narrow Balmer emission profiles without evidence 
for an associated P Cygni wind absorption are probably diffi- 
cult to explain in a spherically symmetric model for the wind of 
HD 108. In fact, the width of the lines indicates that they must 
be formed over a small range in radial velocity and hence the 
radial extent of the line forming region in the expanding wind 
should be rather small. Under such circumstances, a spherically 
symmetric wind would most probably result in a pronounced P 
Cygni absorption component that is not observed. 
On the other hand, the different behaviour of the He II A 4686 
and N III emission lines on the one side and the He I and Balmer 
lines on the other side, clearly indicates that there must be at 
least two distinct regions in the atmosphere of HD 108 where 
emission lines are formed. 

7. Discussion 

The spectrum of HD 108 looks pretty much like that of the 
extreme 08 laf star HD 151804 (Crowther & Bohannan |1997[ ), 
but exhibits additional emission lines (C II, O II,...) and 
presents spectacular H I and He I line variations. 



Table 6. Aspect of the H/3, H7 and He I A 4471 Hnes since 
1919. The codes for the references are the following: P124 
for Plaskett (119241), Me25 for Merrill et al. (119251), Sw42 for 



Underbill (1994) found similarities between the spectrum of 
HD 108 and t he spectra of Ofpe/WN9 stars (Bohannan & 
Walborn 1989| ) and those of B[e] s tars. Comparing with more 



recent data (e.g. Nota et al. 1990, we notice some marked 



Swings & Struve ( |1942D, B e50 for Beals (|1950|), Ma55 for 
Maninno & Humblet (1955), H068 for Houziaux & Ringuelet- 
Kaswalder (\196^ , An73 for Andrillat et al. (|1973| ), Hu75 for 
Hutchings ( [19751 Un94 for Underbill ( 11994| ) and Ba99 for 
Barannikov ( |1999| ). 



Date 


H/3 


H7 


He I 


reference 


1 01 Q 


f Lyg. 




P Lyg. 


iviezj 




F Cyg. 






PI 0/1 


1 m A 1 mo 


P Cyg. 


PCyg. 


Ads. 


riejU 


1 O/i 1 


F Cyg. 


P Cyg. 


ADS. 


iviaj J 


1 QA1 


F Cyg. 


P Cyg. 




OW4Z 










Ma55 


1944-1945 


Abs. 


Abs. 


Abs. 


Be50 


1950 


Abs. 


Abs. 


Abs. 


Ma55 


1953 


Abs. 


Abs. 


Abs. 


Ma55 


1954 


PCyg. 


Abs. 


Abs. 


Ma55 




(weak) 








1966 


PCyg. 






H068 


1968-1974 


PCyg. 


PCyg. 


PCyg. 


An73 










Hu75 


1976 


PCyg. 


PCyg. 


PCyg. 


this work 


1982-1985 


PCyg. 


PCyg. 




Ba99 


1986-1991 


PCyg. 


PCyg. 


PCyg. 


Un94 


1989-1991 


PCyg. 


PCyg. 




Ba99 


1991 




PCyg. 


PCyg. 


this work 


1993-1994 




PCyg. 




this work 


1996-1997 


PCyg. 


PCyg. 


Abs. 


this work 


1998-1999 


PCyg. 


Abs. 


Abs. 


this work 


2000 


Abs. 




Abs. 


this work 



of HD 108 is sharper and the He I lines can evolve into pure 
absorptions, not only weaker P Cygni. 

HD 108 was also discussed in connection to the Luminous 
Blue Variable (LBV) stars. However, until now, HD 108 
has not presented spectacular spectral type variations. The 
He II A 4542 line is only slightly variable and the apparent 
changes in spectral type discussed in Sect. 3.1 rather reflect the 
variations of the wind emission in He I A 447 1 . Finally, unlike 
most LBVs or LBV candidates, no nebular emission was 
detected around HD 108, even if the star has been suggested in 
the past to be the exciting star of a faint nebula (see Hi ggs & 
Ramana 1968 ). But, in a more recent work Lozinskaia ( 1982| ) 
excludes the presence of a nebula around HD 108. 



differences between the spectrum of HD 108 and spectra of 
OfpeAVNL stars. For example, the He II and N III absorptions 
are stronger in the spectrum of HD 108. Also, the Ha emission 



Andrillat et al. (1982) presented some characteristics of Oe 
stars. In their sample of stars. Ha, He I A 6678 and H I Paschen 
lines are often present in emission. These emissions are 
variable, on time scales of several years. They can completely 
disappear or even become absorption lines. HD 108 presents 
similar characteristics of varying He I and H I lines, and its 
near IR spectrum is also similar to those of Oe stars. HD 108 
does not display the Fe I I AA7 515, 7712 and O I A7772 lines 
seen by Andrillat et al. (1982) in most Oe spectra, but these 
lines are also missing in the spectra of some confirmed Oe/Be 
stars, as for instance in X Persei. 
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H? He I A 4471 A 

Fig. 9. Upper panels: aspect of H7, He I A 4471 and Hf3 for some recent years. The lines gradually evolve from P Cygni profiles 
into pure absorption lines. Bottom panels: reconstructed profiles of the same lines (see text). The He II contribution has been 
subtracted from H7 and Hf3 (bottom, left and right) while the He I photospheric component has been subtracted from the He I 
total profile (bottom, middle). The heliocentric correction has not been applied to the spectra displayed in this figure. 



Moreover, Divan et al. (1983) showed that Oe and Be stars 



are brighter and redder when the emissions in the Balmer 
continuum and in Ha are stronger Unfortunately, there are 
few photometric studies of HD 108. Its V magnitude was 
measured by Plaskett & Pearce (1931), Hiltner & Johnson 
( |1956| ), Blanco et al. ( |1968| ) and Leitherer & Wolf ( |1984| ) 
who all found V=7.40 and by Bouigue et al. (1961), who 
found V=7.35. The B-V color index is about 0.18 for all 
articles, except the last one (B-V = 0.14). However, HD 108 
was sometimes attributed a variable character: in the New 
Catalogue of Possible Variable Stars (NVS), its V magnitude 
varies from 7.35 to 7.48. Chilardi et al. ( |1948[ ) reported 
magnitude variations up to 0.4 mag, and more recently. 



Barannikov ( 1999| ) showed that there was an abrupt change in 
the light curve in 1996 and 1997, when He I A 4471 went into 
absorption, followed quite soon by H7. At that time, the star 
became suddenly fainter and bluer (in A(V-R)). In addition, 
the Hipparcos satellite has measured a mean Hp magnitude of 
7.412±0.012 mag with a maximum of 7.39 and a minimum 
of 7.43 (Hipparcos main catalogue, ESA 1997). The small 
amplitude of the light curve argues against HD 108 being a dor- 
mant LBV, but it could be compatible with an Oe/Be behaviour. 

HD 108 seems thus to share several characteristics of Oe 
or Be-type stars. However, the spectrum of HD 108 displays 
rather sharp emissions without any indication of a double- 
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peaked morphology that could be attributed to an equatorial 
disc or flattened wind. Therefore, if HD 108 has indeed an 
equatorial disc, it should be seen under a rather low inclination. 
In fact, spectropolarimetric observations (HaiTies 2000| ) do not 
exclude the possibility of a disc, provided it would be seen 
face-on. 

Moreover, as mentioned in the introduction, several conflicting 
determinations of the mass-loss rate of HD 108 can be found 
in the literature. From the unsaturated resonance lines seen 
in the lUE spectrum of HD 108, Hutchings & von Rudloff 
( 1980) and Howarth & Prinja (1989) inferred rather low values 
of 2 X 10^'' and 3 x 10~^ M q yr~^ respect ively. Using the 
Ha line flux, Peppel ( |l984 and Leitherer ( |l988| ) obtained 
M ~ lO~^M0yr~^. Finally, the largest mass-loss rates 
(5 X 10^^ M0 yr^^ comparable to those of Wolf-Rayet stars) 
were derived from the infrared excess (Ferrari-Toniolo et al. 



1981). The various values listed here refer to observations ob- 



tained at the same epoch (between 1978 and 1982) thus ruling 
out any long-term trend as the origin of the discrepancies. 
These different determinations of M could provide further 
support for a scenario based on a disc seen face-on. In fact, the 
UV resonance lines are sensitive to the absorption along the 
line of sight, whereas the IR continuum probes the emission 
from the extended disc. 

As already quoted in the first section, HD 108 has been 
sometimes classified as a runaway: Bekenstein & Bowers 



(1974) found a peculiar velocity Vp > 98 kms while 



Underhill ( |l994| ) gives Vp — 43 km s ^ . The Hipparcos satel- 
lite has measured the proper motion of HD 108 and found 
Macos((5) = (— 5.39 ± 0.66) X 10~^ arcsecyr"^ and — 
(— 0.93±0.54) X lO^arcsec yr~^. Us ing th e method described 



by Moffat et al. ( 1998 and eiTatum in 1999 ), we calculate : 



(w)pec = (-3.14 ± 1.01) X lO-^arcsecyr-i 
if^b)pec = (0.47 ±0.56) X 10-3— 
yielding {vt)p ec 
from the Sun equal to 2.51 kpc (Gies |1987[ ) and a 30% uncer- 
tainty in the distance. 



arcsec yr 

— 37.8 ± 16.5 kms^^, assuming a distance 



Moffat et al. (1998) define a runaway as a star with {vt)p 



42 



> 

This criterion does not allow us to attribute a 



definite runaway status to HD 108. 
8. Conclusions 

Using a 30 years long campaign of spectroscopic observations, 
we have analysed the strange behaviour of HD 108. The analy- 
sis of our RV time series failed to reveal a significant period 
and our data rule out the periods claimed in previous stud- 
ies (e.g. Hutchings 1975, Barannikov 1999). However, HD 108 
presents some outstanding long-term variations: He I and H I 
lines changed recently from P Cygni profiles to apparently pure 
absorption lines. A similar situation was observed about 56 
years ago. The line profile variations are thus very probably 
recurrent. We tried to remove the underlying photospheric ab- 
sorption to better understand some of the variations. The re- 
stored Balmer and He I lines consist of rather narrow emis- 
sion lines overlying a broad (photospheric) absorption line and 
reveal no evidence for a P Cygni type absorption. The other 



emission lines either display lower amplitude variations (N III 
AA 4634-41, He II A4686,...) or are only marginally variable 
(e.g. the O II emissions, the unidentified AA4486, 4504 emis- 
sions,...). These various emissions have thus most probably dif- 
ferent origins. Several different hypotheses on the wind struc- 
ture and the nature of HD 108 can be considered: 



1. The emission lines in the spectrum of HD 108 could be 
formed in jets as suggested by Underhill (1994). The long- 
term variations could be due to a strong occultation of the 
jets by the stellar body. Such an occultation might result 
from a precession of the jets. In our data, the absorption 
continues to strengthen after the emission apparently van- 
ished (e.g. He I A 4713) whereas the EWs should reach a 
limiting value during the eclipse. One possibiUty could be 
that in 2000 we are still in a pre-eclipse phase. However 
several details remain to be addressed. Where do the rather 
constant emissions (e.g. He II A 4686) originate? What trig- 
gers the precession of the jets? Moreover, the narrow emis- 
sion lines indicate that these jets must be seen almost per- 
pendicular to the line of sight and it seems difficult to imag- 
ine a geometric configuration that could explain the secular 
variations without exhibiting a simultaneous modulation of 
the emission on the (much shorter) time scale of the stellar 
rotation period. 

2. The star could be surrounded by a disc which undergoes 
density oscillations (see e.g. the models by Okazaki 1991 ). 
The emission lines would be attenuated when these oscilla- 
tions are eclipsed by the star, provided that the inclination 
is sufficient. However, the narrow emission lines without 
any evidence of a double-peaked profile and the spectropo- 
larimetric observations of Harries (2000) suggest that such 
a disc should be seen face-on. Therefore it seems difficult 
to explain the observed variations with this model. 

3. Alternatively, the wind of HD 108 could be confined by 
a strong enough magnetic field into a high density cool- 
ing disc in the pl ane of the magnetic equator (see Babel & 
Montmerle 1997). The emission lines would be formed in 
this cooling disc. If, for some reason, the strength of the 
magnetic field changes as a function of time, we expect 
the density in the disc and hence the emission strength to 
vary. Again, we would have to explain why there is a di- 
chotomy in the behaviour of the emission lines and what 
would be the mechanism causing the modulation of the 
magnetic field strength. 

4. HD 108 might be a binary system harboring a compact ob- 
ject which is revolving in a highly eccentric orbit around 
an Oe star which is surrounded by a disc seen face-on. 
Near periastron passage, the X-rays emitted by the accret- 
ing compact object (neutron star or black hole) could dras- 
tically alter the ionization in the circumstellar enveloppe. 
The X-rays could completely ionize hydrogen and He I, 
and the emission in those lines would thus disappear. As 
long as the compact companion is far enough from the Oe 
star, the He II A 4686 formation zone should not be af- 
fected too much. If the compact object follows a highly 
eccentric orbit, the spectral lines should not be affected 
during most part of the orbit. As some older data are 
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missing (especially from 1955 to 1968), we cannot check 
what was the morphology of the lines at this epoch. On 
the other hand, the enhanced interaction around perias- 
tron should generate an enhanced X-ray emission. In this 
context it is worth mentioning that HD 108 was detected 
as a fairly bright source in the ROSAT All Sky Survey 
(RASS). In fact, Rutledge et al. ( |2000[ ) quote a ROSAT- 
PSPC count rate of 0.0614 ± 0.0124ctss-i for HD 108. 
Assuming that the X-ray emission is produced in an op- 
tically thin thermal plasma with kT — 0.5 keV and ac- 
counting for an interstellar neutral hydroge n colu mn den- 
sity of 3.4 X 10^^ cm"2 (Diplas & Savage |l994, we de- 



rive an unabsorbed flux of 2.3 x 10~ ergcm^ s^ in 
the energy range 0.1 - 2.0ke 'Vp|. Again assuming a dis- 
tance of 2.51 kpc (Gies 1987| ), we derive an X-ray lumi- 



nosity of Lx = 1.7 X 10 erg s^ corresponding roughly 
to logLa;/Lboi ~ —6.0. The ratio between the X-ray and 
bolometric luminosities in 1990 (i.e. at the time of the 
RASS) was therefore about six times larger than expected 
from t he 'canonical' relation proposed by Berghofer et al. 
(|1997|). Observations of HD 108 wifli the XMM-Newton 



satellite might allow to confirm the compact companion hy- 
pothesis since they will be obtained when the stars would 
be close to periastron passage. 

If this compact companion model applies, then the sys- 
tem should certainly have undergone a supernova. Since 
the runaway status is not established, another confirmation 
of this scenario could come from the detection of a super- 
nova remnant (SNR) in the vicinity of the star. No such 
SNR has been detected so far. In the ambient circumstellar 
medium ionized by the O-star, the SNR signature would be 
very difficult to see. Another possibility for a null detection 
could be th at this remnant has already evaporated. Maeder 
& Meynet (1994) predict lifetimes between 2.6 x 10^ and 
4.4 X 10^ yrs for stars of spectral type 03 to 06, whereas a 
typical 07.5 star stays 4.6 x 10^ yrs on the main sequence^. 
Since a SNR takes typically 5 x 10* yrs to evaporate, the 
compact companion hypothesis cannot be ruled out on the 
grounds of the lack of a detected SNR. 

As the behaviour of HD 108 seems recuiTent, and could repro- 
duce after about 56 years, this star certainly deserves a long- 
term spectroscopic monitoring. The H I and He I lines could 
either reach a steady state or continue their evolution towards 
stronger absorptions, before they gradually recover a P Cygni 
profile. Photometric observations of this star should also be 
useful to assess the amplitude of the possible magnitude varia- 
tions of HD 108, which might be in correlation with its line pro- 
file variability. Moreover, multiwavelength studies should be 
undertaken, in order to confirm possible variations of the mass- 
loss rate or the compact companion hypothesis. Only then, will 
we finally be able to answer the question about the real nature 
ofHD108. 



This flux is in very good agreement with the 0.032 cts s 



EINSTEIN-IPC count rate quoted by Chlebowski et al. ( |1989[ ). 

' Assuming that the typical initial masses of 03, 06 and 07.5 stars 
are respectively 85, 30 and 20 Mq. 
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